Quantum Non-Magnetic state near Metal-Insulator Transition - a Possible Candidate 

of Spin Liquid State 
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In this paper, based on the formulation of 0(3) non-linear a model, we study two dimensional 
7r — flux Hubbard model at half-filling. A quantum non-magnetic insulator is explored near metal- 
insulator transition that may be a possible candidate of spin liquid state. Such quantum non- 
magnetic insulator on square lattice is not induced by frustrations. Instead, it originates from 
quantum spin fluctuations with relatively small effective spin-moments. In the strong coupling 
limit, our results of spin velocity and spin order parameter agree with results obtained from earlier 
calculations. 
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People have been seeking for quantum spin liquid 
states in spin models with predominantly antifcrromag- 
netic short ranged interactions for over two decades 
For example, various approaches show that quantum spin 
liquids may exist in two-dimensional (2D) S = 1/2 Ji- 
J2 model or Heisenberg model on Kagome lattice. In 
these models, the quantum spin liquids are accessed (in 
principle) by appropriate frustrating interactions. In par- 
ticular, such type of spin liquid states can be described 
by the Hubbard model formalism in the strong coupling 
limit. 

Recent experiments on the triangular lattice show that 
the spin liquid ground state may be realized in the or- 
ganic material k - (BEDT - TTF) 2 CU 2 (CN) 3 @, 0, [|. 
Motivated by experiments, U(l) slave-rotor theory of 
Hubbard model on triangular lattice^ and its SU(2) 
generalization on honeycomb lattice were formulated [6|. 
It is predicted that quantum spin liquids may lie in the 
insulating side of the metal-insulator (MI) transitions. 
Because the spin liquid is adjacent to the MI transitions, 
people may guess it is the local charge fluctuations rather 
than frustrations that disrupt spin ordering and drive the 
ground state to a spin liquid. Such type of spin liquid 
can be described by the Hubbard model formalism of the 
intermediate coupling region. 

Recently, it becomes hot to use ultracold atoms as sim- 
ulators of quantum many-body systemsQ ■ In particular, 
the 7r-flux Hubbard model (or the Hubbard model with <f> 
flux) on square lattice has been designed with ultracold 
atoms in an optical lattice. An artificial magnetic field of 
7r-flux (or <fi flux) in an optical square lattice is proposed 
to be realized by different approaches, such as laser as- 
sisted tunneling method^], laser methods by emplo yin g 
dark states^ or dressing two-photon by laser fields [loj] . 
Without the nesting condition, the MI transition of the 
7r-flux Hubbard model may differ from that of the tradi- 
tional Hubbard model on square lattice. Thus, due to 
nodal fermions in the non- interacting limit, it becomes 
an interesting issue to learn the MI transition of the 
7r-flux Hubbard model. In addition, it is known that the 
insulator state of the 7r-flux Hubbard model belongs to 
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FIG. 1: Illustrations of a Ti-flux lattice. There is a 7r-flux 
phase when an atom hops aroud a plaquette (the gray rect- 
angle). Where a is the length of the side that is chosen to be 
unit. 



a special class of antiferromagnetic (AF) ordered state - 
nodal AF insulator (NAI), an AF order with relativistic 



massive fermionic excitations [11]. Another issue here is 
whether the nodal AF insulator is a long range AF order 
or a short range one. 

In the followings an 0(3) non-linear a model (NLerM) 
is developed to investigate properties of NAI in the 7r-flux 
Hubbard model. Based on the NLcM, we will show that 
a non-magnetic insulator (a short range AF order) may 
exist in the NAI of 2D 7r-flux Hubbard model when the 
spin fluctuations are considered. 

Metal-Insulator transitions of the ir-flux Hubbard 
model - The Hamiltonian of 2D 7r-flux Hubbard model 



is 



H = - 



tijc\cj 



U 



n^mi 



(1) 

Here Cj = (cjf , c«j.) T are defined as electronic annihilation 
operators. U is the on-site Coulomb repulsion, /j is the 
chemical potential and at half- filling it is ^. denotes 
two sites on a nearest-neighbor link, njj and are the 
number operators of electrons at site i with up-spin and 
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FIG. 2: The staggered magnetization of the n-flux Hubbard 
model (circle solid line) and that of the traditional Hubbard 
model (square solid line) at zero temperature. (U/t) cl ~ 3.11 
is the critical point of the metal-insulator(MI) transition of 
the -n-flux Hubbard model. 



down spin, respectively. There is a 7r-flux phase when 
an atom hops around a plaquette in a 7r-fiux lattice(See 
Fig.l). So the nearest neighbor hopping tij in a 7r-fiux 



lattice could be chosen as 12 t 



i,i+y 



te 



±*4 



Because the Hubbard model on bipartite lattices is un- 
stable against antiferromagnetic (AF) instability, at half- 
filling, the ground state may be an insulator with AF or- 
der (NAI). Such AF order is described by the following 
mean field result 



It is well known that due to the nesting effect, there is 
no MI transition of the traditional Hubbard model (arbi- 
trary interaction will lead to a magnetic instability) . The 
situation is difference for the MI transition of 7r-flux Hub- 
bard model. The MI transition of the 7r-flux Hubbard 
model occurs at a critical value about U/t ~ 3.1lfl2l] 
(See Fig.2). In the weakly coupling limit (U/t < 3.11), 
the ground state is a semi-metal (SM) with nodal fermi- 
points[H|. In the strong coupling region (U/t > 3.11), 
due to M 0, the ground state becomes an insulator 
with massive fermionic excitations. By contrast, there 
is only the insulating phase of the traditional Hubbard 
model (See Fig.2). However, the non-zero value of M 
only means the existence of effective spin moments. It 
does not necessarily imply that the ground state of NAI 
is a long range AF order because the direction of the 
spins is chosen to be fixed along z-axis in the mean field 
theory. Thus one needs to examine stability of magnetic 
order against quantum fluctuations of effective spin mo- 
ments based on a formulation by keeping spin rotation 
symmetry. 

Effective nonlinear a model of spin fluctuations - In 
the following parts we will focus on the NAI state and 
don't consider local charge fluctuations and the ampli- 
tude fluctuations of M that are all gapped in the region 
of M ^ o[l3j]. To deal with the spin fluctuations, we 
use the path-integral formulation of electrons with spin 

The interaction 



14, \m m uja m 



rotation symmetry 
term in Eq.([T]) can be handled by using the SU(2) in- 
variant Hubbard-Stratonovich decomposition in the ar- 
bitrary on-site unit vector O,; 



(-l)VM) 



(2) 



Here M is the staggered magnetization. For the cases of 
spin up and spin down, we have a = +1 and a = — 1, 
respectively. Then in the mean field theory, the Hamil- 
tonian of the 2D 7r-fmx Hubbard model is obtained as 



H 



YjtiAtj + h.c.) - J2 ( 3 ) 



where A = leads to the energy gap of electrons and 
o z is the Pauli matrix. After diagonalization, the spec- 
trum of the electrons is obtained as 



A 2 



(4) 



where |£k| = \/4i 2 (cos 2 k x + cos 2 k y ) corresponds to the 
energy of free fermions. By minimizing the free energy at 
temperature T in the Brillouin zone, the self-consistent 
equation of ([3]) is reduced into 



k 



(5) 



where N is the number of the particles. 



-[a 



■cjaci 



(6) 



Here a = (o~ x , cr yi <j z ) are the Pauli matrices. By replac- 



ing the electronic operators c] and Cj by Grassmann vari- 
ables c* and Cj , the effective Lagrangian of the 2D 7r-flux 
Hubbard model at half filling is obtained: 

C e s = 2J c i9 T Ci - ^2 (ti,jC^Cj + h.c.) - A ^2 c*fti-crci. 

i (ij) i 

(7) 

In particular, we describe the vector il^ with the hal- 
dane's mapping: 



n i =(-l)*n i yi-L? + L i . 

Here = (nf ,nf) is the Neel vector with n 2 = 
and Li is the small transverse canting field with L, ■ n; 

00, MM- 



(8) 



Then we rotate to z-axis at each site on both sublat- 
tices by performing the following spin transformation 14, 
El EE E3, 13, *l>i = U iCi , U}xvi ■ aUi = a z and 
Ujhi ■ aUi = U ■ a. After the spin transformation, the 
effective Hamiltonian becomes: 



3 



H cS = 5>;aa (i) ^ - E + ^- c -) -A^^[(-1)V,^1 - If + h ■ <#< 

i 

i 



(9) 



where the auxiliary gauge fields a,-j = ciij. i<j x +ctij,2 and M = ( — 1)* (^*c z ^»j), we obtain the effective Hamilto- 
a (i) — oo,i (i) o-a;+ao,2 («) cr y are defined by e Mij = U-Uj nian: 
and a (i) = lljd T Ui. In terms of the mean field result, 



(10) 



In this equation we have used the approximations 

y/l - If ~ 1 - | and e io « ~ 1 + ia l3 . 

In the next step, we integrate the gapped fermion fields 
and get the quadric terms of [do (?) — Act • lj] and a^. 
Then the effective action becomes: 



1 /' ,3 

S eS = - / dr V[-4<r(a (i)-A<r ■ li) 2 +4/9 s a: 

2 Jo 



2A 2 

2 1 ZAA i2l 

(11) 



lowing two equations 21 



1 E 



n ^2(iar + A2)i 



1 E 



A 2 



^ V 4(i^r + A 2 )^ 



where the parameters p s and ^ are derived from the fol- and the corresponding coefficient e 2 is given as: 



(12) 
(13) 



e 2 = i 2 [cos (2k x ) (A 2 + 8t 2 + At 2 cos (2k y )) + A 2 + M 2 + i 2 cos (4fc x )] 

I 



To learn the properties of the low energy physics, 
we study the continuum theory of the effective action 
in Eg. pip . In the continuum limit, we denote rij, lj, 
idij ~ UjUj - 1 and a (i) = U}d T U l by n(x,y), l(x,y), 
U^d x U (or WdyU) and Wd T U, respectively. From the 
relations between Wd^U and <9 M n, 



2 2,2 
= , + a. 



1 



ao-1 = -- (n x d T n) ■ 1, 



(14) 



Vi + a M:2 = 7(^*1) > ^ = x,y, (16) 



(17) 



the continuum formulation of the action in Eq. ([TT|) turns 



- --(3 r n) , t = 0, (15) into 



1 Z 1 ' 3 /■ 2A 2 
S cS = ±J drj d 2 r[,(d T nj 2 + Ps ( V nf - AtA, (n x d T n) ■!+(—- 4A\)1 2 ] 



(18) 



where the vector a is defined as (00,1, 00,2: 0) • obtain the effective NLctM of the 7r-flux Hubbard model 

Finally we integrate the transverse canting field 1 and 
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FIG. 3: The spin stiffness p a and the spin wave velocity c of 
the 7r-flux Hubbard model. 



Sc« = ^ f dr [ d 2 r[- (d T nj 2 + c( V n) 2 ] 
25 Jo J c 



(19) 



with a constraint n 2 = 1. The coupling constant g and 
spin wave velocity c arc defined as 



9 

,2 



Ps 
Ps 

X ± 



(20) 
(21) 

is the 



Here p s is the spin stiffness and x — ( — 2UY 
transverse spin susceptibility. 

The numerical results of p s and c of the 7r-flux Hubbard 
model are illustrated in Fig. 3, where one can find that 
p s = 0.03936* = 0.2460J, c = 0.226278t = 1.41424J 
in the strong coupling limit (U = 25i) match the earlier 
results p s — -j, c = \/2J (J = tt) that are obtained 

n [H. 



from the Heiscnberg model [22|, [23, 24], 

In addition, we need to determine another important 
parameter - the cutoff A. On the one hand, the effec- 
tive NLcrM is valid within the energy scale of Mott gap, 
2A = UM. On the other hand, the lattice constant is a 
natural cutoff. Thus the cutoff is defined as the following 
equation A = min(l, ^)[l5j]. 

Magnetic properties of the nodal AF insulator : In this 
section we will use the effective NLcrM to study the mag- 
netic properties of the insulator state. The Lagrangian 
of NLcrM with a constraint (n 2 = 1) by a Lagrange mul- 
tiplier A becomes 



off 



1 

2^ 



(c) T rif +c 2 (vn) 2 + U(l-n 2 ) (22) 



where iX = m 2 and m is the mass gap of the spin fluctu- 
ations. 



FIG. 4: The dimensionless coupling constant a of the n- 
flux Hubbard model (cicle solid line) and that of the tra- 
ditional Hubbard model (square solid line). There are three 
regimes, semimetal(SM), quantum disordered (QD), antifer- 
romagnetic(AF), separated by two critical points (U/t) * ~ 



3.11, (U/t) c 



4.26, respectively. There is only the AF 



regime on the traditional square lattice. 



At finite temperature, by rescaling the field n — > V Nn 
and using the large-N approximation, the solution of 
rig = (n) is always zero that is consistent to the Mermin- 
Wigncr theorem. From Eq. (|22j) , we may get the solution 
of m as 



m 



2TsinlT 



sinh( 



cA 
2fc B T 



)]■ 



(23) 



At zero temperature, the solutions of uq and m of 
Eq. (|22p are determined by the dimensionless coupling 
constant a ~ gA. In particular, there exists a critical 
point a c — 47r ( or g c = ^f) : For the case of a < 47r, 
we get solutions of n and m as n = (1 — ^-) : ^ 2 and 
m = 0. For the case of a > 4ir, we get solutions of no 
and m as n = and m = 47rc(^ — i). So we calcu- 
late the dimensionless coupling constant a = gA of the 
7r-flux Hubbard model and show results in Fig. 4. The 
quantum critical point corresponding to a c = An turns 
into U/t ~ 4.26 which devides the NAI state into two 
phases - a quantum disordered state (QD) in the region 
of 3.11 < U/t < 4.26 and a long range AF order in the 
region of U/t > 4.26. The results show sharp contrast 
to those from the traditional Hubbard model, where the 
dimensionless coupling constant is always smaller than 
a c = 47r. 

In the region of U/t > 4.26 (where a < a c ), at low 
temperature the mass gap m of spin fluctuations is de- 
termined by 



27TC 1 1 

m~2fc B Texp[- — ( )]. 

ksT g g c 



(24) 
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FIG. 5: The mass gap m of the spin fluctuations and the 
ordered spin moment Mo of the 7r-fiux Hubbard model at 
zero temperature. There are three regimes, semimetal(SM), 
quantum disordered (QD), antiferromagnetic(AF), separated 
by two critical points (U/t) cl ~ 3.11, (U/t) c2 ~ 4.26, respec- 
tively. 



Because the energy scale of the mass gap m is al- 
ways much smaller than the temperature, i.e., m <C 
k^T (or uj n ), quantum fluctuations become negligible 
in a sufficiently long wavelength and low energy regime 
(m < |cq| < k^T) . Thus in this region one may only con- 
sider the purely static (semiclassical) fluctuations. The 
effective Lagrangian of the NLctM then becomes 



C 



2 



(Vn) 



(25) 



where 



Ps 



(i-i) 



is the renomalized spin stiff- 
ness. At zero temperature, the mass gap vanishes (See 
Fig. 5(a)) which means that long range AF order appears. 
To describe the l ong range AF order, we introduce a spin 
order parameter 26J, [27|, |28| 



M M. g.,,r, 
M = -n Q = - { l-f)^ 



(26) 



As shown in Fig. 5(b), the ground state of long range AF 
ordered phase has a finite spin order parameter, fn addi- 
tion, in the strong coupling limit, U/t — > oo, the values 
naturally match the results derived from the Heiscnbcrg 
model mapped from the 7r-flux Hubbard model. 

In the region of 3.11 < U/t < 4.26 (where a > a c ), 
there is a finite mass gap of spin fluctuations, m = 
47rc(^ — i) at zero temperature (See Fig. 5(a)). There- 
fore, the ground state of the insulator in this region is 
not a long range AF order. Instead, it is a quantum 
disordered state (or non-magnetic insulator state) with 
zero spin order parameter Mo = (See Fig. 5(b)). The 
existence of a non-magnetic insulator state provides an 
alternative candidate for finding spin liquid state. 



FIG. 6: Illustrations of the relations between the coupling 
constant g and the staggered magnetization M of the 7r-flux 
Hubbard model(cicle solid line) and the traditional Hubbard 
model (square solid line in inset (b)). 



Obviously, such type of quantum non-magnetic insu- 
lator on bipartite lattices is induced neither by geome- 
try frustrations regarded as the examples in varied spin 
models nor by the local charge fluctuations with finite 
energy gap. What is the physics origin of this quan- 
tum non-magnetic state? The key point is that, due 
to the special electron dispersion ( the existence of nodal 
fermions for non-interacting case) the coupling constant 
g is almost proportional to jj near the MI transition 
(See Fig. 6). Hence the non-magnetic state originates 
from quantum spin fluctuations of relatively small effec- 
tive spin-moments, M — > 0. 

Let us compare the properties of the insulator state 
in 7T-flux Hubbard model and those in the traditional 
Hubbard model. For the traditional Hubbard model on 
square lattice, due to the nesting effect, there is no MI 
transition at finite U and the insulator state here doesn't 
belong to NAI. In the U/t — ► limit, the coupling con- 
stant g is not proportional to -^(See Fig. 6). Instead, g 
is about g ~ -^(y) 1 / 4 that becomes smaller and turns 
into zero the weakly coupling limit (See more details in 
Ref. [Hj]). So the quantum fluctuations of the effective 
spin-moments are suppressed. Using the NLerM formu- 
lation, due to g < g c (See Fig. 4), the ground state of the 
Hubbard model on square lattice always has a long range 
AF order. 

Conclusion - In this paper, to deal with the spin fluc- 
tuations, we use the path-integral formulation of elec- 
trons with spin rotation symmetry and then the effective 
NLerM is obtained to describe the NAI state of the tt- 
flux Hubbard model. We calculate the spin stiffness, the 
transverse spin susceptibility, the spin wave velocity and 
the coupling constant g. In the strong coupling limit 
(U/t — > oo), our results of spin velocity and spin order 
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parameter agree with the results obtained from earlier 
calculations of the traditional Hubbard model. How- 
ever, we find that the coupling constant g in the NAI 
state of the 7r-flux Hubbard model shows different be- 
haviors to that in the insulator state of the traditional 
Hubbard model. In particular, a quantum non-magnetic 
insulator state (3.11 < U/t < 4.26) is explored near the 
MI transition that corresponds to the strong coupling 
region of the effective NLctM, g > g c . Such type of quan- 
tum non-magnetic insulator in bipartite lattices is driven 
by quantum spin fluctuations of relatively small effective 
spin-moments. 

Such non-magnetic insulator state is differ- 
ent from that proposed in organic material 
K- (BEDT-TTF) 2 CU 2 (CN) 3 by U(l) slave-rotor 
theory in Ref.Q. Firstly, the non-magnetic insulator 
state is really a short range AF insulator (although we 
don't know its exact nature) followed by a long range AF 
order with increasing U; however, the spin liquid state 
in Ref. [H is really a U(l) spin liquid with spinon fcrmi 
surface, of which no long range AF order appears with 
increasing U. Secondly, the local charge fluctuations play 
important role in the slave-rotor theory; in contrast, the 
local charge fluctuations are irrelevant here. Therefore, 
our results illustrate a new candidate for finding spin 
liquid state. 

An interesting issue is the nature of the non-magnetic 
insulator. Is it a valence-band crystal [2j|, or algebra 
spin liquid state 0, [3(|, ...? In addition, another issue is 
whether there exists non-magnetic insulator state of the 
Hubbard model in honeycomb lattice, of which there also 
exist nodal fcrmions. These issues arc beyond the scope 
of the present work and will be left for a future study. 
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